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Abstract 

Exclusive p^p^ production in two-photon collisions between a quasi-real and a 
mid-virtuality photon is studied with data collected at LEP at centre-of-mass en- 
ergies 183 GeV < y/s < 209 GeV with a total integrated luminosity of 684.8 pb~^. 
The cross section of the process 77* p°p° is determined as a function of the pho- 
ton virtuality, Q^, and the two-photon centre-of-mass energy, W^^, in the kinematic 
region: 0.2 GeV^ < Q"^ < 0.85 GeV^ and 1.1 GeV <W^^<?> GeV. 
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1 Introduction 

Recently, the L3 Collaboration measured the processes 77* — >■ p°p° and 77* — > p'^p~, where 
one of the interacting photons, 7, is quasi-real and the other, 7*, is off- mass-shell and has a 
virtuality in the range 1.2 GeV^ < < 30 GeV^ [1)2]. The cross sections of these isospin- 
related reactions have a similar dependence on the two-photon centre-of-mass energy, W^^, 
and are of similar magnitude, though the p^p~ cross section is systematically higher than the 
one. These features of p pair-production at high are in contrast with the observed 
suppression, and different W^^ dependence, of p^p~ production [3] with respect to p°p° [4, 5] 
in the data for ^ and W^^ < 2 GeV. 

The observed behaviour of p pair-production at large momentum transfer is well described 
by the QCD-based model developed in Reference 6, as shown by the analysis of the L3 data 
presented in Reference 7. On the other hand, p pair-production by quasi-real photons is still 
not well understood, despite a wide range of theoretical models [8,9]. Thus, the study of the 

evolution of p pair-production between these two regimes is an important task in the 
experimental investigation of vector meson pair-production in two-photon interactions. This 
Letter presents results on the measurement of the process: 

+ — +— * + — 0( 
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in a kinematic region of intermediate values of the squared momentum transfer: 

0.2 GeV^ <Q^ < 0.85 GeV^ 

and for an invariant mass of the hadronic system, W.yy, in the interval: 

1.1 GeV <W.,^<3 GeV. (3) 

The data sample used was collected by the L3 detector [10] at LEP at centre-of-mass energies 
183 GeV < ^/s < 209 GeV and corresponds to an integrated luminosity of 684.8 pb^^. Scattered 
beam electrons^^ which have radiated photons with virtualities in the range (j2)) can be detected 
("tagged") by the Very Small Angle Tagger (VSAT) [11]. The VSAT is an electromagnetic 
calorimeter made of BGO crystals installed around the beam line on opposite sides of the L3 
detector, at 8.05 m from the interaction point. Its geometrical acceptance covers the polar angle 
range 5 mrad < < 10 mrad, for azimuthal angles in the ranges —1.25 rad <0< 1.25 rad and 
TT— 1.25 rad < 0< 7r-|-1.25 rad. When the electron with the largest scattering angle is detected 
by the VSAT, the maximum virtuality of the two photons, Q^, is, to good approximation, equal 
to the transverse momentum squared, pf, of the final state hadron system: 

Q2 = 2EhEsil - cos Os) ^ E^EsOl ^ , (4) 

where Ei, is the beam energy, and Eg and 9^ are the energy and the scattering angle of the 
tagged electron, respectively. The VSAT provides a means to ensure selection of exclusive final 
states by correlating the scattered electron and the detected hadron system. 

The p^p^ production cross section is determined as a function of PV^^ and Q^. The results 
are compared to the generalised vector dominance model (GVDM) [12]. A measurement of 
process in a similar kinematic region was performed at lower centre-of-mass energy by 
the PLUTO Collaboration [5]. The present measurement represents a tenfold increase of the 
statistics compared to that measurement. 

Throughout this Letter, the term "electron" denotes both electrons and positrons. 
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2 Event Selection 



The reaction ((H), contributing to the process 

e'''e~ — >■ e'''e^J^g7r^7r^7r^7r^ , (5) 

is identified by a scattered electron, Ctag, detected in the VSAT and four charged pions measured 
in the tracking chamber. These events are collected by two independent track-triggers [13]. 
The trigger efficiency is determined from the data itself, making use of the redundancy of the 
triggers, and is around 94%. 

Single-tagged events are selected by requiring one electromagnetic cluster in the VSAT. 
This cluster must have energy greater than 50% of the beam energy in order to reduce the 
background and to ensure a sufficient containment of the electromagnetic shower. 

The event candidates must have exactly four tracks with zero total charge. All tracks must 
come from the interaction vertex, have transverse momentum greater than 100 MeV and an 
energy loss in the tracking chamber compatible with the pion hypothesis. 

Events containing muons are removed from the selected data sample. A search for secondary 
vertices is performed and events with reconstructed neutral kaons are also rejected. Energy 
depositions above 60 MeV in the electromagnetic calorimeter, not associated with a charged 
track, are selected as photons. An event is allowed to contain no more than one such photon 
with energy below 300 MeV which should not exceed 10% of the total energy of the four-pion 
system. Events containing higher-energy photons are discarded. 

According to equation (jU, the transverse momentum squared pi of the four-pion system is 
used to measure the of the event. It is required to be in the range (j21). For an exclusive 
final state, the projections of the momentum vectors of the electron tag and the four-pion 
system on to the plane perpendicular to the beam direction must be back-to-back. Therefore, 
the acoplanarity angle, 0aco, calculated from the difference between the azimuthal angles of 
the tagged electron and the four-pion system, is required to be less than 0.4 rad, as shown in 
Figure [T^. 

After all cuts, 1958 events are observed. Their four-pion mass spectrum is shown in Fig- 
ure Et- The region (jH)) is populated by 1836 events, which are used for the cross section 
determination. The mass distribution of the tt+tt" combinations of the selected events, dis- 
played in Figure Eb, shows a strong signal. A prominent clustering of entries is observed 
at the crossing of the p° mass bands in the correlation plot of the masses of the neutral 7r"'"7r~ 
combinations, shown in Figure Efc- No such resonance structure is observed in the correlation 
plot of the masses of the tt+tt"'" and tt^tt^ combinations, presented in Figure Eli- These features 
of the two-particle mass correlations give evidence for a signal from p^p^ intermediate states. 

We also inspect the two- and three-pion mass distributions in the data for production of 
higher-mass resonances. The only statistically signiffcant evidence is for production of the 
/2(1270) resonance, which appears in the two-pion mass spectra in the intervals 2.1 GeV < 
W^^ < 2.5 GeV and 2.5 GeV < PV^^ < 3 GeV as illustrated in Figure El Measurement of 
/2(1270) production is beyond the scope of the present study, which is concentrated on p^ 
pair-production. 
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3 Monte Carlo Modelling and Studies 



To estimate the number of p^p^ events in the selected four-pion data sample, we consider 
non-interfering contributions from the processes: 



To take into account /2(1270) production in the region W^^ > 2.1 GeV, we also consider 
contributions from the processes: 



Monte Carlo samples of processes (jHl) and (jZj) are generated with the EGPC [14] program. 
About 4 million events are produced for each of the processes ©, about 3 million events for 
the first of the processes (|7j) and 1.6 million events for the two remaining processes. The W^^ 
and dependence are those of the 77 luminosity function [15] and only isotropic production 
and phase-space decays are included. The generated events are passed through the full L3 
detector simulation using the GEANT [16] and GEISHA [17] programs and then processed in 
the same way as the data, reproducing the detector behaviour as monitored in the different 
data-taking periods. The scattered electrons are propagated from the interaction point to the 
VSAT taking into account the influence of the magnetic field of the L3 solenoid and the LEP 
quadrupole magnets installed between the L3 detector and the VSAT [11]. 

For acceptance calculations, Monte Carlo events are assigned a Q^-dependent weight, eval- 
uated using the GVDM form-factor [12] for both photons. The detection efficiencies of process 
(jT)), calculated taking into account the detector acceptance and the efficiency of the selection 
procedure, are in the range of 2% — 4% and are listed in Tables [T] and El in different and W^^ 
intervals. The efficiency is mostly limited by the kinematics of the two-photon reaction which 
boosts the hadronic system along the beam direction, often resulting in low-angle tracks outside 
the fiducial tracking volume. This geometric acceptance is then further reduced by the limited 
angular coverage of the VSAT. The detection efficiencies for the other subprocesses from Q 
are of the same magnitude as the p°p° one and follow a similar evolution with and ly^^. 
Including the /2(1270) branching fraction into two charged pions, the detection efficiencies for 
the 77* — i> /2p° and 77* — >■ /27r+7r~ processes are of the order of 2% and the detection efficiency 
of the 77* 72/2 process is about 1.2%. 

For Monte Carlo events passing the selection, the generated energy of the tagged electron 
always exceeds 90% of the beam energy, with an average < Eg/ Eh > = 0.987. This ensures 
that the approximation of by p^, given by relation (0]), is valid within 1% in the region (j21). 
The resolution is determined by the measurement of pf and varies between 8% and 10%; 
the resolution on W^^ is better than 3%. 

4 Background Estimation 

The contribution to the selected sample from e"'"e~ annihilation is negligible. Using 2 million 
Monte Carlo events of the reaction e+e~ —>■ e+e'r'^r" generated with the program LEP4F [18], 
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the background contribution from this process is estimated to be 0.6 ± 0.3 events and is ne- 
glected. The background is mainly due to partially reconstructed events from two-photon 
interactions with higher particle multiplicities in the final state, when tracks or photons escape 
detection. Another background contribution arises from "fake tags", i.e. random coincidences 
with off-momentum beam electrons, which give a signal in the VSAT. These signals correspond 
to energy depositions comparable with the beam energy, and are thus not removed by the cut 
on the energy of the VSAT cluster. 

To estimate the background due to feed-down from higher-multiplicity final states, we select 
a data sample of doubly- charged four-pion events, 7r+'7r"'"7r"'"'7r~ and 7r+7r^7r^7r~, in which at least 
two charged particles are undetected. In addition, we also select 7r+7r^7r+7r~7r° events which 
are used to account for background events with undetected neutral pions. 

All these events are required to pass the event selection procedure, releasing the charge- 
conservation requirement for the doubly- charged events and considering only the 7r+7r^7r+7r~ 
subsystem of the 7r"'"7r~7r+7r~7r° events. The (paco distributions of the accepted background-like 
data events are combined with the distribution of selected Tr+Tr'Tr+Tr" Monte Carlo events so 
as to reproduce the 0aco distribution observed in data. The result of this procedure, applied for 
the events in the kinematic region defined by Q and (jH)), is shown in Figure ^ The resulting 
background levels are quoted in Tables [T] and El 

Dedicated studies show that the off-momentum beam particles at the VSAT location are 
dominantly on the outer side of the LEP ring. Therefore, the related background would appear 
as an excess in the number of events having a tag on the outer side of the accelerator ring, Ngut, 
with respect to the inner side, A^j„. This feature is observed, for instance, in data when the 
cut on (paco is released, as shown in Figure [Ho. In the selected data, displayed in Figure Hh, the 
ratio Nout/Nin = 1.02 ± 0.05 is close to unity, indicating that the background from fake tags is 
small. The ratio of the number of selected events having tag in the forward versus backward 
directions along the beam line, 1.04 ± 0.05, is also compatible with unity. We note that since 
the two background-like data samples used in the background estimation originate from real 
physics processes, they contain a fraction of events with fake tags and take into account the 
effect of this background. 

5 Fit Method 

In order to determine the differential p°p° production rate, a maximum likelihood fit of the 
data to the sum of the processes © and ((7j) is performed in intervals of and W^^. 

The parameter set, fl, comprising the six two-pion masses in an event, namely the four neu- 
tral combinations tt'^tt' and the two doubly- charged combinations tt^tc^, provides a complete 
kinematic description of a four-pion event in our model of isotropic production and decay. For 
each data event, i, with measured variables f2j, we calculate the probabilities, Pj{fli), that the 
event resulted from the j-th production mechanisms of the six possible ones as listed in (jH)) and 
dZj). A likelihood function is defined as: 

A = nE^.^.(^^)' EA, = 1, (8) 

i j=l j=l 

where the fit parameter Xj is the fraction of process j in the -K^n^n^n^ sample for a given 
or W^^ bin and the product runs over all data events in that bin. The probabilities Pj are 
determined by the six-fold differential cross sections of the corresponding process, using Monte 
Carlo samples and a box method [19]. 



5 



In the fits we assume that the processes involving /2(1270) production contribute only 
for W^^ > 2.1 GeV, as suggested by the spectra in Figure El We find the /2 content in the 
data to be well described by the /27r+7r~ and /2p° contributions only. Therefore, in order to 
reduce the correlations between the fitted parameters and their uncertainties, we exclude the 
process 77* — >■ 72/2 from further consideration. Thus, we perform a five-parameter fit in the 

bins and in the W^^ bins for the region W^^ > 2.1 GeV, whereas the fits in the W^^ bins 
for W.yy < 2.1 GeV have three parameters and take into account only contributions from the 
processes (jHl). 

As a check of the fit method, we find that the maximum-likelihood fit reproduces the 
p^p" content of Monte Carlo test samples within statistical uncertainties. Since the analysis 
procedure is optimised for deriving the p°p° contribution, in the following only the p°p° content 
and the sum of the rest of the contributing processes, denoted as "other in" , are considered. 

To check the quality of the fit, the tt+tt^ mass distributions of the data are compared with 
those of a mixture of Monte Carlo event samples from the processes © and (|7|), in proportions 
determined by the fit. The data and Monte Carlo distributions are in good agreement over 
the entire and W^^ range. As an example tt+tt^ mass distributions are shown in Figure El 
The Monte Carlo production model also provides a good description of the measured angular 
distributions, as shown in Figure jH 

6 Results 

The cross section, AcXge, of the process e"''e~ e~^e~p^p^ is measured in bins of and W^^. 
The results are listed in Tables HI and |2l together with the efficiencies and the background 
fractions. The statistical uncertainties, listed in the Tables ^ and |2l are those of the fit. The 
differential cross section daee/dQ"^, derived from AcTge, is listed in Table (H When evaluating 
the differential cross section, a correction based on the Q^-dependence of the p°p° Monte Carlo 
sample is applied, so as to assign the cross section value to the centre of the corresponding 
bin [20]. 

To evaluate the cross section, a^^, of the process 77* , the integral of the transverse 

photon luminosity function, Ltt^ is computed for each Q"^ and PV^^ bin using the program 
GALUGA [21], which performs 0{a^) QED calculations. The cross section a^^ is derived from 
the measured cross section AcTee using the relation AcTee = LTTC-yy- Thus, a^^ represents an 
effective cross section containing contributions from both transverse and longitudinal photon 
polarisations. The cross section of the process 77* p^p^ is listed in Table [H as a function of 

and in Table |21 as a function of W^^. The sum of the cross sections of the other contributing 
processes is also given in Tables [H and jH 

Several sources of systematic uncertainty are considered. The contribution of the selection 
procedure, as estimated by varying the selection criteria, is in the range 4% — 8%. Monte Carlo 
statistics give a contribution in the range 1.5% — 2.3%. The variations of the acceptance ob- 
served when a p-pole form-factor is used instead of a GVDM form-factor for re-weighting Monte 
Carlo events are in the range 1% — 3% for most of the kinematic region. The uncertainties of 
the trigger efficiency, as determined from the data, are in the range 1.9% — 4%. In order to 
estimate the systematic uncertainty of the fit procedure, the size and the occupancies of the 
boxes in the box-fit are varied, as well as the number of bins in which the data is divided for the 
fits. In particular, the fits in are performed using only three bins, which results in the same 
integrated cross section as in the case of four bins. A contribution of 3% — 7% is derived. 
Finally, an uncertainty of 2% — 4% is associated with the background determination. 
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All contributions are added in quadrature to obtain the systematic uncertainties quoted in 
Tables ID and El 



7 Discussion 

The cross section of the process 77* p^p^ as a function of W^-y is plotted in FigureEK, together 
with the sum of the cross sections of the other contributing processes. The shoulder in the latter 
is due to the contribution of the subprocesses involving /2(1270) production. The measured 
pOpO QYoss section shows a broad enhancement at threshold. Figures |SJd and Efc compare the 
measured cross sections with those measured at high [1]. All cross sections decrease with 
and the variation with of the 77* p^p^ cross section is more rapid for low values of 

The measured differential cross section daee/dQ"^ of the reaction e"^e~ e"^e~p°p° is shown 
in Figure together with the high-Q^ data from Reference 1. It is fitted to a form [22] 
expected from QCD-based calculations [23]: 

dt^ee 1 /q^ 

dQ2 ~ Q"(g2+ < >2)2 ' ^ ) 

where n is a constant and < W^^ > is the average W^^ value of 1.8 GeV for this measurement. 
Although this formula is expected to be valid only for ^ W^^, we find it provides a good 
parametrisation of the evolution of all data in the interval 0.2 GeV^ < < 30 GeV^, with 
an exponent n = 2.9 ± 0.1. In the fit, which results in /d.o.f. = 6.9/10 and is shown by the 
line in Figure IHK, only the statistical uncertainties are considered. 

The measured cross section of the process 77* — > p°p° as a function of is shown in 
Figure IHb, together with the L3 data for p°p° production at high [1] and the PLUTO 
measurement for 1 GeV < W^j < 3.2 GeV [5]. The two data sets agree for > 0.3 GeV^ 
while for low values the L3 data lie below the PLUTO measurement. The L3 data is fitted 
with a form-factor parametrisation based on the GVDM model [12], which is found to reproduce 
well the dependence of our measurements. Only the statistical uncertainties are considered 
in the fit, which results in x^/d.o.f. = 7.5/11. Figure Eb also shows the result of a p-pole 
form-factor fit to the PLUTO data, as in reference 5. The L3 data cannot be described by the 
steeper fall of the p-pole parametrisation. 
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range 




Bg 


A(Tee [ pb ] 


daJdQ' [pb / GeV] 


(j^^ [ nb ] 


(T^^ [ nb ] 


[ GeV^ ] 


[%] 


[%] 








other 47r 


0.20 - 0.28 


2.4 


8 


12.5 ±1.18 ±0.79 


155 ±15 ±10 


9.65 ±0.92 ±0.62 


15.6 ± 1.19 ±0.90 


0.28 - 0.40 


3.7 


9 


10.9 ±0.90 ±0.72 


89.5 ± 7.4 ± 5.9 


8.18 ±0.68 ±0.54 


13.0 ±0.89 ±0.86 


0.40 - 0.55 


3.0 


12 


6.37 ±0.78 ±0.54 


42.1 ± 5.1 ± 3.6 


5.59 ±0.68 ±0.47 


12.7 ±0.94 ±0.90 


0.55 - 0.85 


2.0 


20 


6.80 ±0.95 ±0.83 


22.1 ± 3.1 ± 2.7 


4.63 ± 0.65 ± 0.57 


7.86 ±0.81 ±0.79 



Table 1: Detection efficiencies, e, background fractions, Bg^ and measured production cross sections of the reactions e+e~ — > e"'"e~p°p°, 
77* ~^ s-iid of the sum of the rest of the contributing processes, other 47r, as a function of for 1.1 GeV < W^^ < 3 GeV. 
The values of the differential cross sections are corrected to the centre of each bin. The first uncertainties are statistical, the second 
systematic. 
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6.94 ± 1.08 ±0.77 


8.05 ± 1.25 ±0.89 


7.94 ± 1.43 ±0.86 


1.30 - 1.45 


2.6 
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6.81 ±0.85 ±0.58 


11.8 ± 1.48 ± 1.01 


14.3 ± 1.83 ± 1.28 


1.45 - 1.60 


2.8 
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7.07 ±0.81 ±0.62 


13.5 ± 1.55 ± 1.19 


15.9 ±1.83 ±1.30 


1.60 - 1.75 


3.1 


10 


5.61 ±0.70 ±0.47 


11.8 ± 1.46 ±0.99 


16.4 ±1.77 ±1.24 


1.75 - 1.90 


3.1 


10 


3.56 ±0.57 ±0.36 


8.17 ± 1.32 ±0.83 


18.1 ±1.86 ±1.58 


1.90 - 2.10 


3.1 


11 


3.37 ±0.56 ±0.38 


6.38 ± 1.07 ±0.71 


14.0 ± 1.40 ± 1.13 


2.10 - 2.50 


3.2 


11 


2.25 ±0.44 ±0.27 


2.48 ±0.49 ±0.30 


14.1 ± 1.07 ±0.96 


2.50 - 3.00 


3.1 


11 


0.93 ±0.28 ±0.12 


1.01 ±0.31 ±0.13 


6.85 ±0.70 ±0.61 



Table 2: Detection efficiencies, e, background fractions, Bg^ and measured production cross 
sections of the reactions e"'"e~ — >• e'^e~p^p^, 77* — > and of the sum of the rest of the 
contributing processes, other An, as a function of W^^ for 0.2 GeV^ < < 0.85 GeV^. The 
first uncertainties are statistical, the second systematic. 
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Figure 1: (a) Distribution of the acoplanarity angle, (paco, between the electron and the 
7r"'"7r~7r'^7r~' system for data (points) compared to the four-pion Monte Carlo (open histogram) 
and the background estimated from the data (hatched histogram). The arrow indicates the 
selection cut. The shapes of the Monte Carlo and the background are fixed, and their sum 
is normalised to the total number of events, (b) and (c) Distributions of the azimuthal angle 
of the tagged electron in the selected events, 4>tag, ior electrons in the inner side of the LEP 
ring (in) and, folded over it, distributions for electrons in the outer side of the LEP ring (out). 
In (b) all cuts but the acoplanarity cut are applied and in (c) all cuts are applied and the 
corresponding four-pion Monte Carlo distributions are also shown. 
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Figure 2: Mass distributions for the selected events, (a) Mass of the four-pion system, W^^. (b) 
Mass of 7r+7r~ combinations (four entries per event), (c) Correlation between the lower versus 
higher mass combinations of the tt+tt^ pairs (two entries per event), (d) Correlation between 
the masses of the 7r+7r+ and tt^tt^ pairs. The two-dimensional plots in (c) and (d) have a bin 
width of 50 X 50 MeV^ and the size of the squares is proportional to the number of entries. 
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Figure 3: Mass distributions of vr+vr" combinations (four entries per event) for the three higher 
W^.y intervals and for the total sample for 0.2 GeV^ < < 0.85 GeV^. The points represent 
the data, the hatched area shows the component and the open area shows the sum of the 
rest of the contributing processes. The fraction of the different components are determined by 
the fit and the normalisation is to the total number of events. The plot for the entire Wjj 
range, 1.1 GeV < Wjj < 3 GeV, is a sum of the distributions of all fitted Wjj intervals. 
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Figure 4: Comparison of the data and Monte Carlo angular distributions for the kinematic 
regions (1) and (2): (a) | cos 6p |, the cosine of the polar angle of the p° with respect to the 77* 
axis in the 77* centre-of-mass system; (b) | cos | , the cosine of the polar angle of the pion in 
its parent p° helicity-system; (c) A0, the angle between the decay planes of the two mesons 
in the 77* centre-of-mass system; (d) | cos^^ab |j the cosine of the opening angle between the 
two TT"*" directions of flight, each one defined in its parent p° helicity-system. There are two 
entries per event in (a), (c) and (d) and four entries per event in (b). The points represent 
data, the hatched area shows the component and the open area shows the sum of the rest 
of the contributing processes. The fraction of the different components are determined by the 
fit and the normalisation is to the total number of events. 
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Figure 5: (a) Cross section of the process 77* p^p^ (full points) and the sum of the rest of 
the contributing processes (open points), as a function of W^^ for 0.2 GeV^ < < 0.85 GeV^. 
The bars show the statistical uncertainties. The two sets of points have the same binning and 
some points are horizontally displaced for better legibility. Comparison of the results of this 
measurement and that at high [1] for (b) the 77* —>■ p^p^ process and (c) the sum of the 
rest of the contributing processes. 
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Figure 6: The p^p^ production cross section as a function of Q^: (a) differential cross section 
of the process e+e~ c^c^ p^p^ and (b) cross section of the process 77* p^ p^. The full 
points show the results from this measurement, the open points show the results from the L3 
measurement of p^p^ production at high [1] and the squares in (b) show the results from the 
PLUTO measurement [5]. The bars show the statistical uncertainties. The L3 measurements 
are for the interval 1.1 GeV < W^^ < 3 GcV and the PLUTO measurements for 1 GcV < 
W^-y < 3.2 GeV. The line in (a) represents the result of a fit using the QCD-inspired form of 
equation (9). The solid hue in (b) represents the result of a fit to the L3 data based on the 
GVDM model [12] and the dotted line indicates the result of a fit to the PLUTO data using a 
p-pole form-factor. 
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